The Early Steps of Molecule-to-Material Conversion in Chemical Vapor Deposition (CVD): A Case Study

Transition metal complexes with β-diketonate and diamine ligands are valuable precursors for chemical vapor deposition (CVD) of metal oxide nanomaterials, but the metal-ligand bond dissociation mechanism on the growth surface is not yet clarified in detail. We address this question by density functional theory (DFT) and ab initio molecular dynamics (AIMD) in combination with the Blue Moon (BM) statistical sampling approach. AIMD simulations of the Zn β-diketonate-diamine complex Zn(hfa)₂TMEDA (hfa = 1,1,1,5,5,5-hexafluoro-2,4-pentanedionate; TMEDA = N,N,N′,N′-tetramethylethylenediamine), an amenable precursor for the CVD of ZnO nanosystems, show that rolling diffusion of this precursor at 500 K on a hydroxylated silica slab leads to an octahedral-to-square pyramidal rearrangement of its molecular geometry. The free energy profile of the octahedral-to-square pyramidal conversion indicates that the process barrier (5.8 kcal/mol) is of the order of magnitude of the thermal energy at the operating temperature. The formation of hydrogen bonds with surface hydroxyl groups plays a key role in aiding the dissociation of a Zn-O bond. In the square-pyramidal complex, the Zn center has a free coordination position, which might promote the interaction with incoming reagents on the deposition surface. These results provide a valuable atomistic insight on the molecule-to-material conversion process which, in perspective, might help to tailor by design the first nucleation stages of the target ZnO-based nanostructures.     

Identification of Flavonoids as Putative ROS-1 Kinase Inhibitors Using Pharmacophore Modeling for NSCLC Therapeutics

Non-small cell lung cancer (NSCLC) is a lethal non-immunogenic malignancy and proto-oncogene ROS-1 tyrosine kinase is one of its clinically relevant oncogenic markers. The ROS-1 inhibitor, crizotinib, demonstrated resistance due to the Gly2032Arg mutation. To curtail this resistance, researchers developed lorlatinib against the mutated kinase. In the present study, a receptor-ligand pharmacophore model exploiting the key features of lorlatinib binding with ROS-1 was exploited to identify inhibitors against the wild-type (WT) and the mutant (MT) kinase domain. The developed model was utilized to virtually screen the TimTec flavonoids database and the retrieved drug-like hits were subjected for docking with the WT and MT ROS-1 kinase. A total of 10 flavonoids displayed higher docking scores than lorlatinib. Subsequent molecular dynamics simulations of the acquired flavonoids with WT and MT ROS-1 revealed no steric clashes with the Arg2032 (MT ROS-1). The binding free energy calculations computed via molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) demonstrated one flavonoid (Hit) with better energy than lorlatinib in binding with WT and MT ROS-1. The Hit compound was observed to bind in the ROS-1 selectivity pocket comprised of residues from the β-3 sheet and DFG-motif. The identified Hit from this investigation could act as a potent WT and MT ROS-1 inhibitor.     

Atomistic insights into the selective therapeutic activity of 6-(2,4-difluorophenoxy)-5-((ethylmethyl)pyridine-3-yl)-8-methylpyrrolo[1,2-a]pyrazin-1(2H)-one towards bromodomain-containing proteins

Cross-target effect has been one of the major mechanisms of drug toxicity, this has necessitated the design of inhibitors that are specifically tailored to target particular biomolecules. 6-(2,4-difluorophenoxy)-5-((ethylmethyl)pyridine-3-yl)-8-methylpyrrolo[1,2-a] pyrazin-1(2H)-one (Cpd38) is an inhibitor possessing high inhibition rate and tailored specificity towards bromodomain-containing protein 4 (BRD4). In this research, we used an array of computational techniques to provide insight at the atomistic level the specific targeting of BRD4 by Cpd38 relative to the binding of Cpd38 with E1A binding protein P300 (EP300); another bromodomain-containing protein (BCP). Comparatively, binding of Cpd38 improved the conformational stability and compactness of BRD4 protein when compared to the Cpd38 bound EP300. Also, Cpd38 induced a conformational change in the active site of BRD4 that facilitated a complementary pose between Cpd38 and BRD4 suitable for effective atomistic interactions. Expectedly, thermodynamic calculations revealed that the Cpd38-BRD4 system had higher binding energy (−36.11 Kcal/mol) than the Cpd38-EP300 system with a free binding energy of −15.86 Kcal/mol. Noteworthy is the opposing role Trp81 (acting as hydrogen bond acceptor) and Pro1074 (acting as hydrogen bond donor) found on the WPF and LPF loops respectively play in maintaining Cpd38 stability. Furthermore, the hydrogen bond acceptor/donator ratio was approximately 4:1 in Cpd38-BRD4 system compared with 2:1 in Cpd38-EP300 system. Taken together, atomistic insights and structural perspectives detailed in this report supplements the experimental report supporting the improved selectivity of Cpd38 for BRD4 ahead of other BCPs while providing leeway for the future design of BET selective agents with better pharmacological profile.     

New insights into the anti-hepatoma mechanism of Alisol G-metal ions complexes based on c-myc DNA

This study investigated the anti-hepatoma molecular mechanism of Alisol G, which is an effective component of the Chinese medicine Alisma orientalis, in the presence of metal ions Cu²⁺ and Fe³⁺ based on c-myc DNA. Here, a combination of Alisol G and metal ions (Cu²⁺, Fe³⁺) to augment anti-hepatoma efficiencies of Alisol G has been identified by methyl thiazolyl tetrazolium (MTT) assay. Network pharmacology revealed that c-myc DNA was the potential target of Alisol G with respect to its anti-hepatoma effects. By performing multi-spectroscopic analyses, we showed that the interaction of Alisol G with c-myc DNA was a process of static quenching. The binding constants and thermodynamic constants indicated that a 1:1 complex was formed between Alisol G and c-myc DNA. Moreover, metal ions strengthened the interaction between Alisol G and c-myc DNA. Molecular docking and molecular dynamics simulation further unveiled that the higher binding affinity between Alisol G-Fe³⁺ complex and c-myc DNA as compared to Alisol G-Cu²⁺ complex. This probably resulted from the polarization of metal ions and the structural flexion of Alisol G. The C22-O31-H76 and C18-O32-H77 of Alisol G were key groups in the interaction with c-myc DNA. Addition of metal ion, had greatly changed the c-myc DNA-binding domain of Alisol G while didn’t affect the kinetic stability of the interaction, thus facilitating the insertion of Alisol G into c-myc DNA A-T base pair. Importantly, the DG113 of c-myc DNA was important for its binding to metal ions. Together, our findings suggested that Alisol G in combination with metal ions may be an efficient and promising option for the treatment of liver cancer.     

Conformational Selection Mechanism Provides Structural Insights into the Optimization of APC-Asef Inhibitors

Metastasis is the major cause of death in colorectal cancer and it has been proven that inhibiting an interaction between adenomatous polyposis coli (APC) and Rho guanine nucleotide exchange factor 4 (Asef) efficaciously restrain metastasis. However, current inhibitors cannot achieve a satisfying effect in vivo and need to be optimized. In the present study, we applied molecular dynamics (MD) simulations and extensive analyses to apo and holo APC systems in order to reveal the inhibitor mechanism in detail and provide insights into optimization. MD simulations suggested that apo APC takes on a broad array of conformations and inhibitors stabilize conformation selectively. Representative structures in trajectories show specific APC-ligand interactions, explaining the different binding process. The stability and dynamic properties of systems elucidate the inherent factors of the conformation selection mechanism. Binding free energy analysis quantitatively confirms key interface residues and guide optimization. This study elucidates the conformation selection mechanism in APC-Asef inhibition and provides insights into peptide-based drug design.     

On the Effect of the Various Assumptions and Approximations used in Molecular Simulations on the Properties of Bio-Molecular Systems: Overview and Perspective on Issues

Computer simulations of molecular systems enable structure-energy-function relationships of molecular processes to be described at the sub-atomic, atomic, supra-atomic or supra-molecular level and plays an increasingly important role in chemistry, biology and physics. To interpret the results of such simulations appropriately, the degree of uncertainty and potential errors affecting the calculated properties must be considered. Uncertainty and errors arise from (1) assumptions underlying the molecular model, force field and simulation algorithms, (2) approximations implicit in the interatomic interaction function (force field), or when integrating the equations of motion, (3) the chosen values of the parameters that determine the accuracy of the approximations used, and (4) the nature of the system and the property of interest. In this overview, advantages and shortcomings of assumptions and approximations commonly used when simulating bio-molecular systems are considered. What the developers of bio-molecular force fields and simulation software can do to facilitate and broaden research involving bio-molecular simulations is also discussed.     

Influence of temperature and alloying elements on the threshold displacement energies in concentrated Ni–Fe–Cr alloys

Concentrated solid-solution alloys(CSAs) have demonstrated promising irradiation resistance depending on their compositions. Under irradiation, various defects can be produced. One of the most important parameters characterizing the defect production and the resulting defect number is the threshold displacement energies(E_d). In this work, we report the results of E_dvalues in a series of Ni–Fe–Cr concentrated solid solution alloys through molecular dynamics(MD)simulations. Based on several different empirical potentials, we show that the differences in the E_dvalues and its angular dependence are mainly due to the stiffness of the potential in the intermediate regime. The influences of different alloying elements and temperatures on E_dvalues in different CSAs are further evaluated by calculating the defect production probabilities. Our results suggest a limited influence of alloying elements and temperature on E_dvalues in concentrated alloys. Finally, we discuss the relationship between the primary damage and E_dvalues in different alloys. Overall, this work presents a thorough study on the E_dvalues in concentrated alloys, including the influence of empirical potentials,their angular dependence, temperature dependence, and effects on primary defect production.     

Experimental study on age and gender differences in microscopic movement characteristics of students

Campus security has aroused many concerns from the whole society. Stampede is one of the most frequent and influential accidents in campus. Studies on pedestrian dynamics especially focusing on students are essential for campus security, which are helpful to improve facility design and emergency evacuation strategy. In this paper, primary and middle school students were recruited to participate in the single-file experiments. The microscopic movement characteristics,including walking speed, headway, gait characteristics(step length, step frequency and swaying amplitude) and their relations were investigated. Age and gender differences in the headway-speed diagram and space requirements were analyzed by statistical tests. The results indicated that the impacts of age and gender were significant. There were three stages for the influence of gender on the headway-speed diagram for both age groups. The impacts on students' space requirements were consistent for different age and gender groups. But the impacts of age and gender on free-flow speed were affected by each other. Due to the connection of walking speed and gait characteristics, the comparisons of gait characteristics between different ages and genders were performed to understand the corresponding differences in speed more deeply. The results showed that differences in step length and swaying amplitude between males and females were significant for both age groups. The effect of gender on step frequency was significant for primary students. But for middle school students,whether gender had significant impact on step frequency was not clear here because of the large P-value. Besides, the influence of age on gait characteristics changed with gender.     

水化硅酸钙纳米压痕分子动力学模拟方法优化

针对水化硅酸钙纳米压痕模型忽视了压头与基底之间相互作用的问题，由尺寸差异引起的金刚石压头难以计算的问题，以及Wittmann模型无法得到实际接触面积的问题，提出了新的模型与计算方法．结合分子动力学方法，采用金刚石压头-Wittmann模型基底的组合方式构建无定形态水化硅酸钙纳米压痕试验模型．在建模阶段，考虑到压头模型与基底模型粒子间尺寸差异，提出了等比例替换模型，通过公式推导并就不同尺寸模拟结果验证了等比例替换模型的可行性．在计算阶段，提出了局部前处理的弛豫方法进行模拟．确定最大荷载位置处的接触面积为546 nm2，进而求出水化硅酸钙模型硬度H为0.84 GPa、折合模量Er为30.52 GPa．并通过纳米压痕试验，验证了模拟结果的准确性，证明了模型的科学性，对今后水化硅酸钙(C-S-H)纳米层面的模拟具有重要借鉴意义．     

Structure of Zinc and Nickel Histidine Complexes in Solution Revealed by Molecular Dynamics and Raman Optical Activity

The histidine residue has an exceptional affinity for metals, but solution structure of its complexes are difficult to study. For zinc and nickel complexes, Raman and Raman optical activity (ROA) spectroscopy methods to investigate the link between spectral shapes and the geometry were used. The spectra were recorded and interpreted on the basis of ionic equilibria, molecular dynamics, ab initio molecular dynamics, and density functional theory. For zwitterionic histidine the dominant tautomer was determined by the decomposition of experimental spectra into calculated subspectra. An octahedral structure was found to prevail for the ZnHis₂ complex in solution, in contrast to a tetrahedral arrangement in the crystal phase. The solution geometry of NiHis₂ is more similar to the octahedral structure found by X-ray. The Raman and ROA structural determinations of metal complexes are dependent on extensive computations, but reveal unique information about the studied systems.